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Qh Abstract 

We compute the complete spectrum of the displacement Hessian opera- 
tor Tiifi = —mp™~ 2 V • (/9*V?/>), which is obtained from the confined porous 
medium equation in dimensions N > 2 by linearization along vector fields VV> 
around its stationary attractor p*. Here, p* denotes the Barenblatt profile and 
tp is an element in the weighted Sobolev space iJi that arises in a natural way 
r^ as the tangent plane in the (formal) gradient flow interpretation of the porous 

medium equation for the entropy E{p). On a formal level, the operator % 
is conjugate to the Hessian of E(p) via similarity transformation. We show 
that H can be understood as a self-adjoint operator on H x and find that its 
spectrum is purely discrete. The knowledge of the complete spectrum and the 
explicit information about the corresponding eigenfunctions give new insights 



o 
o 



> 

on the convergence and higher order asymptotics of solutions to the porous 
medium equation towards its attractor p*. More precisely, the inspection of 
the eigenfunctions allows to identify symmetries in R N with flows whose rates 
of convergence are faster than the uniform, translation-governed bound. The 
present work imitates the analogous study of Denzler & McCann [9] for the 
fast-diffusion equation and complements their results in the porous medium 
r> regime. 

X 

1 Introduction 

In this paper, we study the long-time asymptotics of nonnegative solutions to the 
porous medium equation, i.e., 

d t p-A{p rn ) = inR* (1) 

with exponent m > 1. This equation is best known for modeling the flow of an 
isentropic gas through a porous medium; other applications include groundwater 
infiltration, population dynamics, and heat radiation in plasmas (cf. [3TI Ch. 2]). 
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Solutions to (fil) feature different phenomena depending on the degree of the non- 
linearity p m . In the case m — 1, equation (JT1) is the ordinary diffusion (or heat) 
equation. For < m < 1, the diffusion flux mp m ~ l diverges as p vanishes and thus, 
for suitable initial data, the solution spreads over the whole R^ immediately. In 
this situation, equation ([I| is often called the fast diffusion equation. In the porous- 
medium range m > 1, the diffusion flux increases with the density and degenerates 
where p = 0. Consequently, solutions will preserve a compact support and hence 
this type of propagation goes by the name slow diffusion. In the present paper, we 
restrict our attention exclusively to the latter case. For a study of the fast diffusion 
equation and more general evolutions of porous-medium type, we refer to Vazquez 
301 and references therein. 



The Cauchy problem for the porous medium equation is solved in various settings. 
As the equation is degenerate parabolic, that is, it is parabolic only where the so- 
lution is positive, solutions are in general not classical. More precisely, if the initial 
datum is zero in some open subset of H N , then there is a slowly propagating free 
boundary that separates the region where the solution is positive from the region 
where it is zero. For suitable initial configurations, e.g. < po £ L 1 ($l N ), unique 
strong solutions are known to exists, and these solutions are bounded and contin- 
uous. Moreover, strong solutions preserve total mass, \\p(t, ■ )\\l 1 (r n ) — WpoWl 1 ^) 
for all t > 0. The Cauchy problem for the porous medium equation is reviewed by 
Vazquez in 151] . 



It is well-known [35] |3], El] that the porous medium equation allows for self-similar 
solutions, so-called Barenblatt solutions, propagating on the scale \x\ = t a where 



a 



iV(m-I) + 2' 



and given by 
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P .(«,x) = j BS L-^-^_) (2) 



Here, (•) + = max{ • , 0} and L is an arbitrary constant that can be fixed, for instance, 
by normalizing the total mass or by choosing the radius of the support of p*. It is 
interesting to note that although the Barenblatt solution is a strong solution of the 
porous medium equation, it is not a solution of the Cauchy problem as p* diverges 
to the delta measure So (times a constant depending on L) at time zero — a reason 
for which it is often called a "source-type solution" . 

The Barenblatt solution describes the long-time behavior of any solution p having 
same total mass as p*. Indeed, for arbitrary non- negative initial data in L 1 , solutions 
spread with a finite propagation speed over the whole space and, while diffusing, the 
shape of the solution becomes asymptotically close to the profile of the Barenblatt 
solution, 

p(t, ■ ) w p*(£, • ) as t > 1. 



This long-time behavior was intensively studied over many years, starting with the 
work by K amin [T5l H~6] , who established uniform convergence in one space dimen- 
sion. The generalization to several dimensions is due to Friedman & Kamin [12] 
and Kamin & Vazquez [IT]. These authors use similarity (rescaling) transforma- 
tions and identify the limit orbits for the porous medium equation via compactness 
arguments. A completely different approach relies on energy (or entropy) methods, 
a way prepared by Ralston in [27], who established a first convergence result in L 1 
in one dimension. This approach was generalized independently by Carrillo & 
Toscani, Otto, and Del Pino & Dolbeault in P [221 [7] to arbitrary dimen- 
sions. In these articles, instead of using geometric properties of the density function, 
the authors study the abstract energy landscape and compute the decay rates of the 
energy (or entropy) functional. In [6j and [22], these decay rates are converted into 
bounds on the asymptotics in L l and in the Wasserstein distance, respectively. We 
refer to Vazquez' excellent review article [29] for an almost complete discussion of 
the long-time asymptotics of solutions to the porous medium equation, including a 
discussion on the optimal rates of convergence. 

The investigation of the optimal rate at which solutions converge to the Barenblatt 
solution is tied together with the choice of the initial data. The minimal requirement 
of integrability of the initial data — which is from the physical point of view natural 
when thinking of po as a concentration density - - is mathematically necessary as 
solutions to non-integrable initial data, i.e., p *= L\ oc \L l , show a different behavior 
in the long-time limit (cf. [291 Cor. 1.2]). The authors in [HI [221 [7] derive algebraic 
rates of convergence for finite entropy solutions, cf. ^ below. The exponent in 
the algebraic rate is estimated by a, which is optimal for spatial translations of the 
Barenblatt solution. The goal of the present article is to go one step further and to 
study the complete spectrum of a linearized version of the porous medium equation. 
The knowledge of the spectrum and of all eigenf unctions not only yields information 
about the slowest converging modes (like a spectral gap estimate would do), but 
also allows to extract information about the characteristic geometric pathologies 
to all orders. In particular, we will study the role of translations, shears, and 
dilations for the convergence rates of solutions towards its attractor p*. While 
spectral properties of the heat equation (i.e., (pj) with m = 1) are well-know in any 
space dimension, the investigation of the spectrum of the porous medium equation 
is completely new in the multi-dimension case. In the one-dimensional setting, 
Angenent [2] computed the long-time asymptotics to all orders for solutions with 
compactly supported initial data. An earlier formal spectral analysis is due to 
Zel'dovic & Barenblatt [34] . Progress in the analogous problem for the fast 
diffusion equation was recently obtained by Denzler, Koch & McCann [8]: the 
authors compute higher order asymptotics in weighted Holder spaces, building up 
on earlier results for the linearized fast diffusion equation derived by Denzler & 
McCann [9]. 

Our approach to compute the spectrum of the linearized porous medium equation 
follows closely the work of Denzler & McCann [9], who studied the analogous prob- 



lem for the fast diffusion equation. Their work is based on a method which is fairly 
common in the quantum mechanics literature, more precisely in the spectral analysis 
of Schrodinger operators, see e.g. [21]. The main difference between the analysis in 
[9] and the one in the present work relies on the occurrence of a free boundary in 
the porous medium regime. In particular, opposed to the situation for m < 1, the 
Barenblatt solution has compact support in 1R, N , This fact will play a crucial role in 
the definition of the linearized operator as ( "asymptotic" ) boundary conditions on 
the boundary of the Barenblatt solution have to be taken into consideration. Some- 
how surprisingly (and unexpected, see [9j p. 303]), compared to the analog study 
in [9] for the fast diffusion regime, the analysis of our linearized operator simplifies 
thanks to the compact support of p* in many aspects. 

We caution the reader that our linearization of the porous medium equation per- 



formed in subsection |1.3| is only formally justified. In particular, the rigorous results 
obtained for the linearized operator have to be considered as conjectures for the 
asymptotics of solutions to the nonlinear equation ([I]). 

1.1 Rescaling 

In order to investigate the asymptotic behavior of solutions to the porous medium 
equation, it is convenient to rescale variables in such a way that the Barenblatt 
solution becomes a stationary solution, that means |x| ~ t a in view of (pi). One 
convenient choice of new variables is 

x = /3t a x, t = exp{a~H), and p = (a/3 2 ) 1/(m-1) *~ aJV P, 



with 



, 2mL 
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v a(m — 1) 
to the effect of 

Pfo X ) = ^N P [ a ln *> Jp 

We remark that the logarithmic rescaling of the time variable, t = ami, guarantees 
the parabolic structure of the equation. In fact, with its new variables, equation (01) 
turns into the nonlinear Fokker-Planck equation 

dip-V-(xp)-A(p m ) = 0, (3) 

sometimes also called the confined porous medium equation. Moreover, with that 
particular choice of (3, the rescaled, time-independent Barenblatt solution p* con- 
centrates on a ball of radius one around the origin, 

' m uvr 1 = (i-m + - (4) 



m 



This normalization will simplify the notation in the subsequent sections. In the 
sequel, we will often refer to p* as the Barenblatt profile. 



The advantage of working with (J3J) instead of (IT]) relies on the fact that the above 
rescaling comes along with a change of perspective: the Barenblatt solution now 
becomes a fixed point of the equation. In fact, as we will see below, for fixed total 
mass M, it is both the unique stationary solution of ^ and the ground state of the 
physical entropy 

E(p) = p(x) m dx+- / \x\ 2 p(x)dx. (5) 

Tfl l J L J 

In this regard, the porous medium equation ([I]) shows a substantially different behav- 
ior than the confined porous medium equation (pi): There is no stationary solution of 
([I]) under the mass constraint ||p||z,i = M for any M > 0. In other words, solutions 
to the porous medium equation do not relax in finite time and the ground state can 
be attained only asymptotically in the long-time limit. 

Somehow related to this feature is the fact that, opposed to the original equation 
(OJ), the confined porous medium equation (J3l) is no longer invariant under spatial 
translations. More precisely, in the new variables the origin is the unique attraction 
point in R w around which the mass asymptotically concentrates, driven by the 
convection term —V • (xp). Due to this loss of translation invariance it is therefore 
not surprising that those initial data po that differ from the ground state p* only 
by a shift, po(x) = p*(x — x ) for some x G 1R, N , will play a distinguished role 
in the discussion of our result on refined asymptotics of ^: Spatial translations 
correspond to the smallest eigenvalue of the linearized operator and generate thus 
those solutions that saturate the optimal rate of convergence. 

The above rescaling, however, also has a small drawback: It is singular at time zero. 
To compare "initial data", one would better consider the moment t — 1. More 
precisely, up to spatial dilations both in the x and p variables, the initial datum 
pit = 0, ■ ) corresponds to the original solution at time one p(t = 1, ■ ). Because of 
our interest in the long-time dynamics, this skew correspondence, however, has no 
influence neither on the subsequent analysis nor on the interpretation of our result. 

It remains to discuss the stationarity of p* for the nonlinear Fokker-Planck equation 
PJ and its minimality for the entropy (J5J) . Calculating the energy dissipation along 
trajectories of the flow ([3]), 



d - 



dt 



xr + 



dx, 



2 m — 
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we see that the dissipation rate is zero if and only if ^\x\ + ^z±p m is constant, 
that is, if and only if p is of the form Q, provided that p and p* have the same 
total mass. This implies that p* is the unique stationary solution of pi). To deduce 
the minimality of p* for the entropy E, we additionally observe that the entropy is 
strictly convex on the convex configuration space of nonnegative densities with fixed 
total mass M, and admits thus a unique minimizer. 



1.2 Otto's gradient flow interpretation and the Wasserstein 
distance 

In his seminal paper [22], Otto introduced a new and physically meaningful gra- 
dient flow interpretation of the porous medium equation (fTl). This interpretation 
is based on the Lyapunov approach of Newman & Ralston in [201 ETJ. In the 
abstract setting, a gradient flow is the dynamical system on a Riemannian manifold 
(Ai,g) that describes the evolution as the steepest descent in an energy landscape, 
i.e., 

d t p + gmd E(p) = 0. (6) 

Here, the gradient grad E of the energy E is defined via Riesz' representation theo- 
rem 

diSE(p).5p = g p (gradE(p),5p) for all p e T p M, 

where di&E denotes the differential of E and T p Ai is the tangent plane at the point 
p G M.. A short formal computation shows that ^ can be equivalently stated as 

d t p = argmin < -g(5p,5p) + diSE(p).5p : 5p G T P M 

In the following, we identify the gradient flow ingredients Ai, g, and E that consti- 
tute the porous medium equation using Otto's approach. (In fact, different gradient 
flow interpretations were proposed for the porous medium equation. Otto's gradient 
flow is natural in view of the thermodynamical background of the equation. More- 
over, it also pertains to the confined evolution, whereas the traditional approaches 
do not.) 

As solutions to (QJ) preserve non- negativity (in fact, solutions satisfy a comparison 
principle) and the total mass M, the manifold in the gradient flow interpretation of 
the porous medium equation has to be chosen as 

M = I p : p > and pdx = M\ . 

Tangent fields Sp at p G A\ must be mass-preserving and non-negative where p = 0. 
However, as the evolution (IT]) does not allow for nucleation of mass outside of the 
support of p, i.e., supp(<9 4 p) C supp(p), it is natural to think of the tangent plane as 

T p Ai = < 5p : supp(<5p) C supp(p) and / 5p = 

Observe that this restriction of tangent fields does not suppress the motion of the 
free boundary since Sp may be nonzero on <9supp(p). The metric tensor is defined 
by 

g p {Sp,5p) = p\Vip\ 2 dx, 
6 



where ip and 5p are related via the elliptic boundary value problem 

-V • (pVip) = 6p in supp(p), _ 

pVip-v = on <9supp(p). 

Here, v denotes the outer normal vector on <9suppp. Actually, arguing even more 
formally, Otto set the elliptic equation in the entire space 1R N , cf. [221 eq. (8)]. 
However, in order to motivate our later analysis, here, we equip the equation with 
its natural boundary conditions. That is, assuming enough regularity of p and 8p, 
the distributional solution ip of —V ■ (pVip) = 5p satisfies ([7|). Again, we remark 
that Vip may be non-zero on <9supp(p), so that the boundary is indeed "free" fo 
move. Equivalently, we can characterize the metric tensor variationally: 

■ygpi&Pi&P) = su p\- p|V<p| 2 cfe - / <p5pdx 

cf. [23j eq. (2.9)], where the supremum is taken over all smooth functions tp on R w . 
The optimal cp is this formulation satisfies ([7]). By identifying tangent fields 5p with 
these new variables, the tangent plane at p can be rewritten as the weighted Sobolev 
space 

T p M = lip : f p\Vip\ 2 dx < oo\ . (8) 

In particular, the metric tensor is finite on tangent fields. 

The energy functional in the gradient flow interpretation is given by the entropy 

E(p) = -±—Jp m dx. (9) 

A short computation now shows the porous medium equation (flj is indeed the 
gradient flow for (M.,g) and E defined as above, see [221 Sec. 1.3] for details. 

As for the confined porous medium equation pj), we remark that p still evolves 
according to the gradient flow on the same Riemannian manifold as p, now with 
total mass M and denoted by (Ai,g), but for the energy E defined in ^. 

We finally address the induced geodesic distance on the Riemannian manifold (Ai, g). 
Benamou & Brenier [5] discovered the relation between the Kantorovich mass 
transfer problem and continuum mechanics by identifying the geodesic distance on 
(Ai, g) between two configurations po, pi having finite second moments J \x\ 2 pi dx < 
oo with the Wasserstein distance 



d>2(pi,Po) = inf // l^i-xol p,(dxidx ), 

H€T(p u p )JJ 

see also [22] Sec. 4.3] or [231 Prop. 4.3]. Here the infimum is taken over the space 
r(pi,po) °f an nonnegative measures p on R^ x R w having marginals p\ and po, 
i.e., 

((xi)p(dxidxo) = / C(xi) Pi(xi)dXi, 



for all ( G C c (R d ) and i = 0,1. For a detailed study of Wasserstein distances 
and more general optimal transportation distances, we refer the interested reader to 
Villani's two monographs [321133] . 

1.3 Linearization 

Instead of linearizing the confined porous medium equation (J3J) around the Baren- 
blatt profile p*, we imitate the approach of Denzler & McCann [§] and compute 
the Hessian of the entropy E at p*. In view of the gradient flow formulation described 
in the previous subsection, this is formally equivalent to linearizing (J3| directly, but 
it has the advantage that a natural class of perturbations is intrinsically given by 
the tangent fields. 

Tangent fields Sp are mass-preserving and concentrate on the support of p*. Fol- 
lowing Otto, and as described in the previous subsection, we identify such tangent 
fields with new variables ip via the elliptic problem (J7|), i.e., 

-V ■ (p*VV>) = Sp in Si, 

p^Vip-u = ondB u { ' 

where B\ denotes the ball of radius one centered at the origin, and thus B\ = 
supp(p*). Moreover, ip is such that 



/ 



p*\Vip\ 2 dx < oo, (11) 



cf. (J8j). Geodesic curves {p s }\s\<^i on M. passing through p* and pointing in the 
direction 5p are obtained by McC Ann's [Tj5] displacement interpolation 

p*(x) = det (i + sV 2 ip(x)) p s (x + sV^(i)] (12) 

(see also [22j Sec. 4.3]), that is, p* is pushed forward by the flux id + sVip( ■ ). We 
note that 

$P = 9 s \ s=0 p s = -V ■ (p*Vtp). 

Recalling that the Hessian hessi? of a function E on a. Riemannian manifold (Ai,g) 
can be calculated as the second derivative of E along geodesies, we define and have 



Hess-E(p„)(V>,^) := g^(Sp,hessE(p^)8p) - 2 



d 2 , . 

s=0 



A formal computation yields: 
HessE(p*)(^,^) 



m 



/pr a («.UP.) a ^+/(^pr 1 +^r)fl?Up.^- 



Observe that supp d 2 \ s=0 p s C suppp* and J d 2 \ s=0 p s dx = 0, and so the second 
integral on the right vanishes thanks to the definition of the Barenblatt profile p*. 
Because of d s \ s=0 p s = — V • (pYV^J, we can rewrite the Hessian as 

Hess£(p*)(^) = m f p™' 2 (v • (p*V^)) 2 dx. 

Finally, integration by parts yields 

Ress E(p,){ilj,ip) = / p*Vip ■ V (-mp™- 2 V ■ (p*W>)) dx. 

At this point, we consider the derivation of the Hessian on a purely formal level. 
The above computations certainly hold true for functions ip that are smooth up to 
the boundary In this case, the integrals are well defined and we can integrate by 
parts since p™ _1 vanishes on <9supp(p*). 

We finally complete this subsection by introducing the displacement Hessian 

-Hij = -m P ™- 2 V ■ (p*Vij) . (13) 



In view of the explicit formula Q for p*, the displacement Hessian can be rewritten 
as 

TTL — 1 -s -s 

nip(x) = (l - |x| 2 ) Aip(x) +x- Vip(x) for x G B x . (14) 

With the help of H, the linearized confined porous medium equation reads dfjj + 
Hip = 0, or on the level of 5p: 

didp + C^UCdp = 0, (15) 

where the operator £ is defined via C5p = ip and 5p and ip are related in the usual 
way. Moreover, a formal analysis yields that % and hess are similar in the sense 
that hess = C~ 1 %C 

The remainder of the paper is exclusively devoted to the analysis of the displacement 
Hessian. We will show in Proposition [l] below, that it defines a self-adjoint operator 
on a suitable domain V(l-L). Our main result, stated in Theorem nl below, gives the 
complete spectrum of H. 



2 Rigorous part 

From now on we claim mathematical rigor. As we focus exclusively on the confined 
porous medium equation (J3J) we may simplify our notation and drop the hats from 
the rescaled variables. Moreover, we restrict to the case 

N>2. 



Let Hp t denote the class of all locally integrable functions on Bi = suppp* such 



that (11) holds, i.e., 

\ 2 H i '■= / p*\Vip\ 2 dx < oo, 

with the identification of two functions if they only differ by a constant. Observe 
that this space is a separable Hilbert space for the topology induced by the norm 
|| • ||#i . In the appendix, we study some properties related to elliptic theory in H\. 

So far, the derivation of the Hessian of the energy functional at the Barenblatt profile 
and the definition of the displacement Hessian were only formally justified. Hence, 
as a first step before debarking into spectral analysis, we have to build a rigorous 
basis for our investigations. Since the formal computations of the preceding section 
are certainly valid for functions that are smooth up to the boundary of the support 
of the Barenblatt profile, i.e., C°°(B~i) functions, let us denote the displacement 



Hessian defined in (13) by W\ Caa i^y A short computation yields that this operator 
is nonnegative and symmetric (cf. proof of Proposition [I]) . In Subsection 2.2 below, 
we show that 'H\c<^ib ) extends naturally to a nonnegative self-adjoint operator % 
with domain 

V{U) = {^E HUB,) n Hi : 

pr -1 V ■ (p*W>) G Hi, p^ ■ v = on dB x } . (16) 

Here, the condition p*Vip ■ v = on dB\ means that 

/ p* V£ • V^ dx = - / f V • (p*W) dx for all f G H l pt . (17) 

We also remark that the regularity that is assumed in the definition of T>(T-C) implies 
that Hip = —mp™' 1 ^ ■ (p^Vip) for all ip G T>(H), where the derivatives have to be 
understood in the weak Sobolev sense. 



Remark 1. The boundary conditions in (17) are asymptotic boundary conditions 
in the sense that 

lim pJxWiptx) ■ -^-r = 0, (18) 

Mti \x\ 

for every function ip G C°°(Bi). In fact, thanks to the density of C°°(B~i) functions 
in Hp t by LemmaMin the appendix, and the Sobolev embedding H^ C L 2 [B\, p 2 ~ m dx) 
in Lemma [3] in the appendix, it is enough to restrict condition (17) to functions 



£ G C°°{B\). Then, the equivalence of (17) and (18) follows immediately (for in- 
stance, by an indirect argument). 

2.1 Main results 

We are now on the position to state our main results: 

10 



Theorem 1. The operator ~H : V(7i) —¥ H^ is self-adjoint. Its spectrum consists 
only of eigenvalues, given by 

N 
X ek = £ + 2k + 2k(k + £+--l)(m-l), 

for £, k G N U {0} with (£,k) 7^ (0,0). The corresponding eig en functions form an 
orthogonal basis of H 1 and are given by polynomials of the form 

1 N iV / x \ 

ipemk(x) = F(-k, + £+ — -l + k;£+—; \x\ 2 )Y £m — \x\ e , 

m — 1 2 2 \\ X \J 

m G {1, . . . , Mi}, where F(a, b; c; z) is a hypergeometric function and Yi m is a spher- 
ical harmonic, corresponding to the eigenvalue £(£ + N — 2) of — A s n-i with multi- 

piidty Mi = {N+e -^r 2) . 

The eigenfunctions corresponding to the eigenvalues A^ are polynomials of degree 
£ + 2k, which are harmonic and homogeneous if k — 0. In fact, the hypergeometric 
functions F(a, b; c; z) in the statement of the theorem reduce to polynomials of degree 



k in the case a = —k, see Subsection 2.4 below. Spherical harmonics Yt m are 



harmonic homogeneous polynomials of degree £. These functions are ubiquitous in 
mathematical physics, see e.g. [13J to name but a single reference. 

We observe that the eigenvalues of the displacement Hessian H, are affine functions 
of the parameter m. There are countably many constant levels A«) = £■ Above 
Aio and A 2 o, the spectrum features a crossing of eigenvalues when varying m. A 
first level crossing occurs for the eigenvalues A30 = 3 and A01 = 2 + N{m — 1) at 
N{m— 1) = 1, and the number of crossings increases with £. The eigenvalues coincide 
with the eigenvalues found by Denzler & McCann [H] for the fast diffusion regime 
in the sense that each eigenvalue continues to the range m < 1 if m is sufficiently 
close to 1. For smaller m, the eigenvalues dissolve into continuous spectrum (that 
disappears in the limit m|l). The occurrence of a continuous spectrum for the fast 
diffusion equation is analytically related to the fact that in the fast-diffusion regime, 
the number of moments possessed by the Barenblatt profile is finite --in contrast 
to the situation for the porous medium equation. In the limit m \. 1, we recover 
the spectrum of the Ornstein-Uhlenbeck operator —A + x ■ V. This coincidence is 
not surprising since the Ornstein-Uhlenbeck operator is the limit operator of ~H as 
mil. In this case, p* is a Gaussian. 

The knowledge of the spectrum of the displacement Hessian % suggests the asymp- 
totic expansion to leading order of solutions p = p(t) to the confined porous medium 
equation close to the stationary solution p*. Arguing purely formally again, we re- 
call the operator C defined by C5p = ip where 5p and tp are related in the usual 



way (10), cf. Subsection 1.3, and observe that the spectrum of V, coincides with the 
spectrum of C~ 1 'HC. Moreover, the corresponding eigenfunctions 5pik m = C^ x ^i.km, 
form an orthogonal basis of the Hilbert space L 2 {B\, p™~ 2 dx). For solutions p{t) 
close to p*, equation ^ reads to leading order 

d t (p(t)-p*) « -C- l UC{p{t)-p*), 

11 
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Figure 1: 2D plots of the supports of the affine transformations (dashed and dotted 
lines) of the Barenblatt profile (solid lines) generated by the eigenfunctions corre- 
sponding to the eigenvalues Aio, A20, A30, and A01 (from the left to the right). 



cf. (15), and thus, exploiting the knowledge of the complete spectrum of T-L (and 



thus of C~ l UC) 
p(t) - p* 

= 2J C £km,5p 



lkm& 



-^tkt 



tkni 



£km i'k'm! 



Ctkml'k'm' 



5p, 



'iJft*m-e" (AftWw)t + 



where the c's are constants that depend on the initial data p(0) only. In order to 
investigate the precise asymptotics of the nonlinear equation (|3]), one has to take 
into account terms of higher order in p(t) — p*. This problem will be addressed in 
future research. In this context we should also mention the work of Koch |TH1 Ch. 
5.3.4], which offers a framework which seems to be suitable for a rigorous treatment 
of higher order asymptotics. For the one-dimensional porous medium equation, the 
full asymptotics was rigorously computed by Angenent in [2j, based on an earlier 
investigation by Zel'dovic & Barenblatt [31]. This present work may be con- 
sidered as a first attempt to complement Angenent's result in the multivariable case. 
Recently, Denzler, Koch & McCann studied the higher order asymptotic be- 
havior of solutions to the fast-diffusion equation |8], using subtle dynamical systems 
arguments in weighted Holder spaces. 

The main purpose of understanding the full asymptotic expansion of solutions in 
terms of the eigenfunctions of the Hessian around the self-similar solution is to get 
a deeper understanding of the long-time behavior of solutions. Although we have 
to keep in mind that our investigation is not fully rigorous in the derivation of 
the displacement Hessian operator, we will in the following discuss the the role of 
the smallest eigenvalues and their corresponding eigenfunctions for the convergence 
towards the attractor p*. We recall that the geodesic curves in the Wasserstein 
space passing through the Barenblatt solution p* along the vector field Vip are 
characterized by p*(x) = det (I + sV 2 ip) p s (x + sVip) (cf. (12)). With this formula 



at hand, the interpretation of the eigenfunctions as generators of mass transport 
is immediate. The smallest eigenvalue Aio = 1 corresponds to transitions along 
the axis e m with m G {1, . . . , Mi = N} (see Figure [IJi). Indeed, in this case, the 
eigenfunction is a homogeneous polynomial of degree one, and the displacement is 
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generated by affine perturbations of the identity map x + sc m e m (here, c m denotes 
a normalizing factor). This fact is in good agreement with the optimal convergence 
rate of solutions of the porous medium equation to the Barenblatt profile, 

d 2 {p{t),p*) < e _t d 2 (po,P*), 

here stated in the Wasserstein topology [221 Theorem 1]. In fact, spatially translated 
Barenblatt profiles saturate this bound. The eigenfunctions to the second smallest 
eigenvalue A20 = 2 generate affine transformations on 1R, N modulo rotations (Figure 
fib). To see this, we observe that ip20m is a harmonic homogeneous polynomial of 
degree two, and thus x + sVip = x + sA 2m x for some symmetric and trace-free 
matrix A 2m . In other words, solutions to (J3J) generated by affine transformations 
show a rate of convergence of precise order e~ 2t to the Barenblatt profile. Denzler 
& McCann [TU] give an exact description of the invariant manifold corresponding 
to these affine perturbations. The role of the next term in the expansion depends 
on the value of (m — 1)JV. As indicated above, there is a first level crossing of the 
eigenvalues A^ and A i: It is A i < A^ precisely for £ > 1/a. Eigenfunctions to 
the eigenvalue A i corresponds to isotropic dilations of the Barenblatt profile (see 
Figure [l]i). Indeed, the spherical harmonic Yqi degenerates to a constant and ip n 
is proportional to 1 — (a(m — l)N)~ l \x\ 2 . It follows that x + sVip = (1 — cs)x for 
some constant c. Such solutions converge with a rate e~ l l a . The mass transport 
corresponding to the (£, 0)-modes with £ > 3 is less obvious, as the corresponding 
eigenfunctions if>^ m are homogenous harmonic polynomials of degree £. We only 
discuss the situation where £ = 3. In this case, the transformation is quadratic 
and leads to triangular deformations (see Figure fit). The rate of convergence is 
e~ 3t . The deformations corresponding to the remaining eigenvalues can be similarly 
computed although the complexity of the underlying symmetries is increasing. 

We finally remark that the knowledge of spectral properties of the displacement 
Hessian 7i, and of the smallest eigenvalue A10 = 1 in particular, immediately yields 
the sharp spectral gap estimate 

p*\Vif)\ 2 dx < Hess E(p* )(■?/>, if)), 

for all if) G Hp t . This estimate, was already derived by Otto [221 Sec. 4.4] and 
builds on McCann's displacement convexity [19] . 



We prove Theoremfljin the remainder of this article: In subsection 2.2 we show that 



the displacement Hessian % can be understood as a self-adjoint operator on H^ 

(Proposition [I]). Moreover, we show that its spectrum is purely discrete provided 

that % has infinitely many eigenvalues that accumulate at infinity (Proposition!^. 



The proof of the latter fact requires some preparation. In Subsection 2.3 we make 
a change of variables into spherical coordinates, which has the advantage that the 
eigenvalue problem reduces to a one-dimensional problem for the radial components 



of the displacement Hessian. This eigenvalue problem is solved in Subsection 2.4 We 
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find that the radial components of H have indeed countably many eigenvalues (and 
so has H). These eigenvalues and the corresponding eigenfunctions are explicitly 
computed in Proposition [3j 

2.2 Self-adjointness and spectral decomposition of 7i 

Since H\ Cao iB-i) is nonnegative and symmetric (cf. proof of Proposition [l| and be- 
cause C°°(Bi) is densely contained in H^ by Lemma p\ in the appendix, the operator 
'H\c<=a(B ) * s closable in H x and its closure is still nonnegative and symmetric. We 
denote this closure again by H. In the following we show that H is self-adjoint with 



domain V(H) defined in (16). 



The following proposition contains a first statement on the spectrum of H, namely 
that the spectrum is real and contained in the ray (0, oo). The first fact is a con- 
sequence of self-adjointness, the second one of nonnegativity together with the fact 
that H is one-to-one and onto (and thus must be in the resolvent set). 

Proposition 1 (Self-adjointness). The operator H '■ T>(H) — > H 1 is nonnegative, 
one-to-one and onto. In particular, H is self-adjoint. 

In fact, it can be shown that H is the only self-adjoint extension of H\c°°ib{)i c ^- P> 
Lemma 4]. For the theory of unbounded operators in general and self-adjointness in 
particular we refer to [281 Ch. 13]. 

The following proposition anticipates some of the results in Proposition [3] below, 
namely that H has countably many eigenvalues that accumulate at infinity. Sup- 
posing this, we combine some properties of H with general results from functional 
analysis to deduce that the spectrum of H is purely discrete. As a consequence, it is 
enough to solve the eigenvalue problem in order to compute the complete spectrum. 
This will be our concern in the two subsections that follow. 

Proposition 2 (Discrete spectrum). Suppose that the displacement Hessian H : 
V(H) — >■ Hp t has countably many eigenvalues that accumulate at infinity. Then the 
eigenfunctions of H form an orthogonal basis of H l . In particular, the spectrum of 
H is purely discrete, i.e., it consists only of eigenvalues with finite multiplicity and 
which do not have a finite accumulation point. 

We start with the 

Proof of Proposition^ Since C°°(B~i) is dense in H l by Lemma p\ and C°°(B\) C 
V{H) C Hp t , we immediately have that H is densely defined. We first show that H 
is a nonnegative, symmetric operator on T>(H). Indeed, by the definition of H and 
by integration by parts we have 

J p^Vcp ■ V {Hil>) dx = m J p™- 2 V • (P* Vp) V • (p*V^>) dx 
for all cp,ij) G V(H). Notice that the boundary term vanishes when integrating by 



parts because of (17) and Hip G H\ . Now nonnegativity follows from choosing 



p* 
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(p — ip and symmetry comes from the symmetry in p and ip of the right hand side 
of the above identity 

Notice also that a densely defined, symmetric operator is automatically one-to-one 
and self-adjoint if it is onto, cf. [281 Theorem 13.11]. In our situation, the latter is 
a consequence of Lemmas M & vol in the appendix. Indeed, given u G H^ , we may 
without loss of generality assume that f pl~ m u dx = as we identify functions that 
are the same up to a constant. By Lemma M we have that u G L 2 (Bi, p1~ m dx). We 
define u = —p 2 ~ ra u and observe that this function satisfies the hypotheses of Lemma 
p5| Consequently, there exists a function ip G Hf oc n if* such that 



-V ■ (p*V^) = u in B u 
p^Vij) ■ v = on dB x . 

By the definition of u, the equation in the ball can be rewritten as — mp™~ 2 V ■ 
(p*VV>) = u, i.e., %?/> = it. Higher interior regularity estimates (cf. [281 Theorem 
8.12]) moreover show that ip G H^ oc fl H^. Hence ip £ T>(H). This proves that "% is 
onto, which completes the proof of Proposition [l] D 

We finally prove Proposition [2] 

Proof of Proposition [5| We first show that the first statement implies the second 
one, and thus we assume that the eigenfunctions of % form an orthogonal basis 
of Hp t . Then the operator % can be written as 7i = J2 n ^n'Pn, where the A n 's are 
the eigenvalues of H and the P n 's the orthogonal projections onto the corresponding 
eigenspaces. We define dE = ^ n V n S n with S n being the delta measure concentrated 
at the point A n , that is, for every Borel set A C R, we have 



E{A) = [ dE(X) = V V n . 



\ n eA 
One easily checks that E defines a spectral measure. Moreover, since 



/+oo 
XdE(X), 
-oo 



the spectral measure E is the unique spectral decomposition of % (cf. [2H1 Theorem 
13.30]). Now the discreteness of E implies the discreteness of the spectrum of "H. 

To prove that H^ admits indeed an orthogonal basis of eigenfunctions of H, it is 
enough to show that % has a compact self-adjoint inverse H~ l : H^ — > H l . In 
fact, in view of the spectral theorem for compact operators on separable Hilbert 
spaces, cf. [2S1 Theorem 4.25], the spectrum of H,~ x is at most countable. Since the 
eigenvalue problems for Hr x and H, are equivalent with reciprocal eigenvalues and 1-L 
has countably many eigenvalues with accumulation point at infinity by assumption, 
the spectrum of Ti^ 1 must also be countable, and thus, the eigenfunctions form an 
orthogonal basis, cf. [2SI Theorem 12.29]. 
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We now show that % has a compact self-adjoint inverse. As % is one-to-one and 
onto, its inverse H _1 : H* — > H^ is well-defined, that is, given £ G if^, we define 
'H _1 £ = ip, where ip is the solution of 



L- 



-m P r 2 V-(p,V^) = £ infii, 
p^V'?/' • ^ = on 9Si. 

In the definition of the Hilbert space H x , we identified functions that only differ by 
a constant. In the following, we fix constants by additionally requiring that 

pl~ m ipdx = ! p 2 ~ m Zdx = 0. 



Notice also that H 1 is bounded and self-adjoint because % is self-adjoint, one-to- 
one and onto by Proposition [IJ cf. [28, Theorem 13.11]. 

It remains to show that T-L" 1 is a compact operator. The argument is standard, 
and is based on a Rellich lemma, namely Lemma [4] in the appendix. We display 
the argument since we are dealing with weighted (and thus non-standard) Sobolev 
norms. Let {£ n }neN denote a bounded sequence in H*. Since H' 1 is a bounded 
operator, also the sequence {?/> n }neN with H' 1 ^ = ip n is bounded. Hence, because 
Hp t is a separable Hilbert space, there exist £ and ip in H* and subsequences (that 
we do not relabel) such that {£ n }neN and {^ n }neN converge weekly to £ and ip, 
respectively. By uniqueness of the weak limit, it is % -1 £ = ip. Moreover, by Lemma 
4hn the appendix, this convergence is strong in 1?{B\, p 2 ~ m dx). Integrating by parts 
and using the definition of ip n , we have 



/ p*\Vlpn\ 2 dx = — / p\ m ip n i n dx. 



By the strong convergence in L [B\,p\ m dx) we conclude that 

lim- I ' p 2 - m ip n i n dx = - I pl~ m iPidx = [ p*\ViP\ 2 dx, 
nfoo m J m J J 

and the last identity follows from integrating by parts again and using H~ l l; = 
ip. This shows that the H^ t norms of the sequence {ip n }ivm converge, and to- 
gether with weak convergence, this yields convergence of {ip n }ntm ( or equivalently 
of {H~ 1 C,n}ne¥i) in Hp t . We deduce that 'H~ 1 is compact. □ 

2.3 Separation of variables in spherical coordinates 

In this subsection, we most widely summarize Subsection 4.1 of [9J. As our argumen- 
tation for the porous medium equation follows closely the one for the fast- diffusion 
equation by DENZLER & McCANN, we entirely skip proofs and refer to [9j for 
details. An important step in the work of Denzler and McCann is the change of 
perspective in the spectral analysis of the displacement Hessian 1-i which comes 
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along with the transformation of the operator into spherical coordinates. In fact, 
the change of variables into spherical coordinates is motivated by the crucial in- 
sight that the displacement Hessian % and the spherical Laplacian A s n-i (i.e., the 
Laplace-Beltrami operator on the sphere S^ -1 ) commute. This is basically a con- 
sequence of the fact that the evolution commutes with rotations. Consequently, 
both operators can be simultaneously diagonalized. In particular, as the spectrum 
of A $ jv-i is well-known, our spectral analysis will reduce to the study of the radial 
part of 7i, which amounts to a one-dimensional operator. We remark that this 
Ansatz is frequently used in the spectral analysis of Schrodinger operators, cf. [21] . 

We transform x G R w into spherical coordinates, i.e., x = ru> with (r, u) G [0, oo) x 
S^ -1 , and recall that under this transformation, the Laplacian reads 

N-l„ 1 



iiRV = d r -\ d r + —A S N-i. 



Here A$jv-i denotes the spherical Laplacian. Since the Barenblatt profile p* is a 
radially symmetric function, i.e., p* = p*(r), one readily checks that T-L and A s jv-i 



commute, cf. (14) 



The eigenvalues of the spherical Laplacian A s n-i are \ii = £{£ + N — 2) with multi- 
plicity Mi = + WJV-2M ~ where M = 1. That is, 

-A S N-iY im = ntYim for £ G N U {0} and m G {1, . . . , M e }, 



and the eigenfunctions Yn m are the spherical harmonics Ye m : S N x — > R, which 
form a complete orthonormal basis for L 2 (S 7V_1 , du): 

Y im (u)Y eml {u)du} = 5 u >5 mm i. (19) 

-i 

For more details about the spectrum and the eigenfunctions of the spherical Lapla- 
cian, see also [131 Ch. 3]. 

Exploiting the orthonormality of the spherical harmonics, we observe that for every 
radially symmetric function / : (0, 1) — > K, the quantity H/^mllH 1 is independent 
of the particular choice of m G {1, . . . , Mi}, 

\m m \\ 2 Hk = f Q ({f\r)Y + ^{f{r)Y)p*{r)r N -'dr. 

This leads us to the definition of the norm ||/||#i = H/y^mll// 1 and the corresponding 
Sobolev space 

H\ = {/ G HUO, 1) = ||/||h]<oo}. 

In the case £ = (i.e., fii = 0), we have to identify functions that only differ by 
a constant. In fact, H\ is a Hilbert space and there is an isometric Hilbert space 
isomorphism 

oo M e 
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given by 



oo M e 
1=0 m=\ 

fim(r) = / ifj(ru))Y em (u})du}, 
where the series are converging in L 2 (S N ~ l , du), and the isometry reads 

oo M e 

\w\h = EE"/«x*- 

1=0 m=l 

The orthogonal decomposition of the Hilbert space H* into eigenspaces generated by 
the spherical harmonics permits us to expand the displacement Hessian into a series 
of radially symmetric operators. More precisely, choosing ip G D(l-C) with i/j(roj) = 
f( r )Yem(w), we have l-Lip = (Hef)Ye m , where %i is the orthogonal projection of H 
onto Hj, namely 

{Htf) (r) = -mpT' 1 (f"(r) + ^V(r) - ^f(r)\ + rf(r). (20) 

It is clear that this operator is again non-negative and self-adjoint with domain 
V(Hi) ={f E H\ : fY em G V(H)}. For an explicit definition of T>{H t ), w e have 
to investigate how the asymptotic boundary condition p*Vip ■ v — behaves under 
the transformation into spherical coordinates. (This investigation is not necessary 
in the work of Denzler and McCann as they consider functions that are spread all 
over K, . In particular, in the fast-diffusion case the Barenblatt profile is positive 
everywhere.) The integrability condition on the free boundary dB\ will become a 
selection criterion for identifying eigenfunctions of Tie, and therefore, it has to be 
studied with care. For ip(rco) = f(r)Yg m (u)) and £ G H l pt) we first compute that 

- / V-(p*V^) dx = - / (p*fr N ~ 1 )' g £m dr + p £ / pJg im r N ~ 3 dr, 
J Jo Jo 

/ p*V^ • V£dx = / p* (f'gffyn + ^fgem) r N ' 1 dr, 
where gimij") = Ljv-i £(ru>)Ye m (u)) dco. Consequently, (17) becomes 

/ p,fg'r N - x dr = - I (pJ'r N ^)'gdr for all g G H\. 
Jo Jo 

As the boundary term at r = automatically vanishes for every f,g G Hj, we simply 
write p*f\ r=1 = to indicate that (21) holds. Now the domain of self-adjointness 
of He can be written as 

V(Ht) = {feH? oc (0,l)nH}\HtfeHl P J'\ r=1 = 0}. 



and 



(21) 



Remark 2. In analogy to the statement (18) in RemarkUi we like to point out that 
the asymptotic boundary condition in the radially symmetric variables reduces to 

limp*(r)/'(r) = 0, (22) 

rfl 

for every function f G C°° [0, 1] . 

2.4 The eigenvalue problem for %£ 

Under the assumption that % has countably many eigenvalues accumulating at 
infinity, we have shown in Proposition [2] above that the spectrum of % is purely 
discrete and that the eigenfunctions form an orthonormal basis of H\- According 
to the isomorphy established in the previous subsection, it is enough to remove this 
assumption for the radial components of %. In the following, we solve the eigenvalue 
problem for the self-adjoint operator Up. 

Proposition 3 (Radial eigenvalue problem). For every leN, the eigenvalue prob- 
lem 

U e f = Xf 

in T>(Hf) has exactly the eigenvalues 

N 
\ k = £ + 2k + 2k(k + £+--l)(m-l), 

where k e N such that (£,k) ^ (0,0). The corresponding eigenfunctions are the 
polynomials 

fa(r) = r £ F(-k } — ^ + £ + j - 1 + k; £ + j; r 2 ), 

where F(-k,b;c;z) = 1 + En=i ( ~(tj"i? n -^ mth 0)« = s(s + I) . . . (s + n - 1). 

From the study of the operator % in Proposition [T] we already know that any eigen- 
value of Hi must be real and positive, i.e., A > 0. Notice that in view of the explicit 
formulas (20) & Q for He and p*, the equation 1-Lif = Xf reads 



,„ , ( N-l 2 r \ ( 2A 1 /iA 

/ + 7-, o )f + \ H 2 ~ ~2 M = m ' X ' 23 

\ r m — 11 — r z J \m — ll — r z r z J 

This is a linear differential equation of second order with singularities at the end- 
points of the interval (0, 1). It is of Fuchsian type, that means, all singular points, 
here 0, 1, and oo, are regular. For every given A, this differential equation has 
two linearly independent families of solutions. It is well known that the study of 
second-order Fuchsian ODEs with three regular singular points is intertwined with 
the study of hypergeometric functions F(a, b; c; z) defined by 

n=0 {C)n 
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where a,b,c,z G R and c is not a non-positive integer. Here, the notation involves 
the Pochhammer symbols (or extended factorials) 

(s) n = s(s + 1) . . . (s + n — 1), for n > 1, and (s)o = 1. 

It is easily verified that the series converges if \z\ < 1. In the case that moreover 
c > b > 0, the hypergeometric function F(a, 6; c; z) has the integral representation 



F(a, b; c; z) 



rfc) 



r(6)r(c-6) 



t b -\l + t) a - c (l + t - zt)- a dt, 



were T(s) denotes Euler's gamma function. We finally quote the fact that is most 
relevant for our purposes: For every choice of a, b,c,z G R, c not a non-positive 
integer, the function F(a, b; c; z) satisfies the Hypergeometric differential equation 



F" + 



a + b — c + 1 
\-z 



F' - 



ah 



zl-z) 



in (0,1), 



(25) 



where the primes indicate derivatives with respect to z, i.e, F' = d z F and F" = d\F . 
For detailed discussions of hypergeometric functions, we refer to [251 ESI H]- A 
compact catalogue of the most important facts can be found in p]. 



Now, we come back to the differential equation (23) and its relation to the hyper- 
geometric equation (25). Our goal is to transform (23) into the hypergeometric 
differential equation (25), for which all solutions and their asymptotics at the sin- 
gular points 0, 1, and oo are well-known. Making the Ansatz f(r) = r 1 F{a ) b; c; r 2 ) 
and supposing that / is a solution, (23) transforms into 

'27 



F" + 



+ 



+ iVl 1 1 

2 z m — 1 1 — z 
7(7 + iV - 2) - m 1 



F' 

7' 



A 



z 2 2(m-l)z(l-z) 



(26) 



with z 



Comparing this differential equation with (25) and identifying a, b, 



c, and d yields a first solution to (23). A second linearly independent solution can 
be deduced from the well-known (but complex) theory of hypergeometric functions. 
We have the following 



Lemma 1. Let 








1/1 „ N \ 1 r- 

1 l P l 11 Ti 


(27) 




2\m-\ U ' 2 l ) 2 vg ' 


b = 


1/1 „ N \ 1 r - 

1 1 P 1 111 ffi 


(28) 


2U-1 2 l ) ' 2 V<? ' 


Q = 


2A / iV \ 2 iV-2 1 

T+ \ i+ ^- 1 ) + T + 7 T^' 

m — 1 \ 2 J m — 1 (to — 1J 


(29) 


c = 


<+£ 


(30) 
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Then the differential equation (23) has two linearly independent solutions fi(r) and 



/bO")- A first solution is of the form 

fi(r) = r e F(a,b;c;r 2 ), 



where F denotes the hyper geometric function defined in (24). In the case where c is 
not a positive integer, a second solution f 2 is given by 



Mr) 



2-l-N 



F{a + 1 - c,b + 1 - c;2 - c;r 2 ). 



If c is a positive integer, c = k G N ; then a second solution has the asymptotics 

fz(r) ~ 
and 



In r for k = 1 . 

r X+I Ik f or k>2\ 



f 2 ( r ) ~ r e+l - 2k as < r < 1. 



As the second solution f 2 will be discarded as a potential eigenfunction of the eigen- 



value problem (23) in the proof of Proposition 3 below, it suffices at this point to 



quote only its asymptotic behavior at zero for positive integers c. 



Proof of Lemma^ A first solution to (23) can be obtained by identifying a set of 
values 7, a, b, and c for which (26) turns into the form (25), and setting /i(r) = 
r 1 F{a 1 b; c; r 2 ). The form of a second linearly independent solution depends on the 
particular values of a, b, and c. In the simplest case, the second solution is obtained 
by an appropriate change of the dependent and independent variables. In some 
cases, however, the construction of the second solution relies deeply on the theory 
of hypergeometric functions, and we have to refer to the relevant literature. 

We start identifying values for 7, a, b, and c such that the hypergeometric function 



F(a,b;c;z) solves (26). Comparing (26) and (25), we immediately see that 7(7 + 
N — 2) = fi£. In view of the definition \xt = £(£ + N — 2) this enforces either 7 = £ 
or 7 = 2 — N — £. We choose 7 = £. Moreover, the constant c is determined by 
7(= £) and iV only, namely by c - 
can be derived from the identities a + b 



N 



- y, i.e., (30). Finally, the values for a and b 



c+1 



m—l 



and ab 



e-\ 



2(m-l)' 



A short 



computation yields the two solutions a± = b T = \ (^zi + ^ + y — l) ± \y/q-> where 



2A 



T + (^+f-l) 2 + 



+ 



1 



JV-2 

m-1 " r (m-1) 2 ' 

= b + and b = a + 



and this quantity is positive since A is 
b-, i.e., p7h-p9b. Observe that c is 



q 

positive. We choose a 

not a non-positive integer since N > 2 and £ > 0, so that F(a, b; c; z) is well defined 

for a, b, and c as above. 

Depending on the particular values of a, b and c, we will in the following compute a 



second linearly independent solution of (23). This solution will again be of the form 
/ 2 (r) = r e F(a,b;c;r 2 ), where F := F(a,b;c;r 2 ) is a second linearly independent 
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solution of the hypergeometric equation (25). Indeed, computing the Wronskian 
W(f h f 2 ) yields ' 

w(f u f 2 ) = Ah-hft 

= 2r 2e+1 (F'F-FF') 
= 2r 2l+1 W{F 1 F) ^ 

since F and F are linearly independent. 

In our derivation of F, we follow [U Ch. 8]. If c is not a positive integer, then 
a second linearly independent solution to the hypergeometric differential equation 



(gjjp is given by z 1 ~ c F(a + 1 - c, b + 1 - c; 2 - c; z) (cf. [2 eq. (8.2.6)]) and thus 
/ 2 (r) = r 2_£ ~ 7V F(a + l — c, 6+1 — c; 2 — c; r 2 ) is a second linearly independent solution 



to (23). Notice that /i coincides with / 2 in the case c = 1 and the latter function 
is not even defined for larger integer values of c. For integers c = k > 1, a second 
solution to the hypergeometric differential equation is given by 

G(a,b;k;z) (log z) F (a, b;k;z) 



F(a + 1 - k)T(b + 1 - k)(k - 1)1 

OO / x ,j\ 



(a) B (6) n 0(a + n) + 0(6 + n) - 0(n + 1) - 0(fc + n) n 



n=0 



(k) n n\ r(a+l-k)T(b + l-k)(k-l)\ 



(_!)*(* - 2)! ^ (a + 1 - fc)„(6 + 1 - k\ k 






r(a)r(6) f- (2 - k) n n\ 



n=0 



r'( s ) 



where #(s) = -p^y and the convention that the last term is zero if k = 1, cf. [U eq. 
(8.4.4)]. We remark that _F(a,6;/c;0) = 1. Inspecting the asymptotic behavior at 
the origin and eventually redefining G(a, 6; k; z) by multiplying by T(a) and/or F(b) 
if a and/or b is a nonpositive integer, we see that 

r-it r r \ f lnz for fc = 1 i 

G(a,b;k;z) ~ j ^_ A for ^ > 2 J as < * < 1. 

and 

d z G(a,b;k;z) ~ z ~ k as < 2 < 1. 

(See also the discussion on page 275 in [4].) It remains to set fi(r) = r e G(a, b; k; r 2 ) 
and notice that k = 1 precisely if £ = and N = 2. This proves Lemma [TJ □ 

We are now in the position to present the 

Proof of Proposition^ We consider the eigenvalue problem T-Lif = Xf for / e 
T>(Hz). By the discussion in the previous section and Proposition [TJ we already 
know that every eigenvalue A must be positive. It is readily checked that the eigen- 
value equation for He is equivalent to the differential equation ( J23J ) . This equation 
is solvable for every A > by Lemma [TJ and a solution /to (23) is an eigenfunction 
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of Tig if and only if / G T>(%^). Moreover, / can be written as a linear combination 



of the linearly independent solutions f\ and f 2 - As solutions of (23) are automati- 
cally smooth in (0, 1) by interior regularity results, we have to study the asymptotic 
behavior of f\ and f 2 at the singular points and 1 in order to decide whether 
fiJ 2 eV(n e ). 

For further references we quote that 

ctb 
d z F(a,b;c;z) = — F(a + 1, 6 + 1; c+ 1;A 
c 

cf. [U eq. (8.2.3)], and therefore 

f x (r) = lr l - x F{a, b; c;r 2 ) + 2— r e+1 F(a + 1,6 + 1; c+l;r 2 ), (31) 

and, provided c is not a positive integer, 

f 2 (r) = (2 -£- Ny- e - N F{a-c+ 1,6- c+ 1,2 -c;r 2 ) 

+ 2 (a-c+l)(b-c+l) r3 _ e _ NF{a _ c + ^ 6 _ c + 2 . 3 _ c . r 2 )(32) 

Moreover, one easily computes that 

F(a,6;c;0) = 1 (33) 

if c is not a non-positive integer. 

Step 1. It holds f 2 G" H\, and thus f 2 is not an eigenfunction ofHu. 

We first consider the case where c = k is an integer. Then, by Lemma [T] it is 

f 2 (r) ~ r ^+!- 2fc asymptotically at the origin, and thus, using 2k = 2£ + N by ( |3~0| 

^(^(/^W) 2 ^" 1 ~ r 2^1-4fe+iV = r l-2* as0<r «l. 

The term on the right is not integrable for any fc > 1, and thus H/^IIh 1 = oo. This 
shows that f 2 (jL H\ in the case where c = k is an integer. If c is not an integer, then 



our argument is based on (32) and (33). In fact, denoting by C a ,b,c a constant that 
depends only on a, b, and c, we have f^r) « (2— £— N)r 1 ~ e ~ N +C a ^, c r 3 ~ e ~ N ~ r x ~ l ~ N 
asymptotically at the origin. Here we have used that £ + iV ^ 2 because c ^ 1. 
Therefore 

P*{r){f2{r)fr N - 1 ~ r 1 ^-^ = r 1 "^ as < r < 1. 

Again, the right hand side is not integrable since c > 1 and thus f 2 $. H\ in the case 
where c is not an integer. 

Step 2. For —a G" N U {0} ; it holds f\ ^ V(7ie), and thus f\ is not an eigenfunction 
ofU e . 

To rule out solutions /i as eigenfunctions of % in the case —a £ NU {0}, we investi- 
gate the asymptotics at the singularity r = 1. Since solutions /i are smooth in (0, 1), 
the asymptotic boundary conditions in the definition of V(7ie) can be explicitly 
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checked via (22). In fact, eigenfunctions must necessarily satisfy lim r |i p*{r)f'(r) = 
0. 

In the following, we show that this condition fails if —a ^NU {0}. In the case that 
a + b — c $. Ti, we have the linear transformation formula 

/ , n T(c)r(c — a — b) , , 

F(a,b;c;z) = ) -^F{a,b;a + b - c+ 1; 1 - z) 

I (c — a)l [c — b) 

+ {l-zr a - b V{ ^^ h ~ C) F{c-a ) c-b ] c-a-b + l ] l-z) ) 



cf. [H 15.3.6]. In view of Q and a + b - c £ Z, it is 

t-,/ 7 \ r(c)r(c — a — b) ,„ .. „ ft r(c)r(a + 6 — c) 

r(c — ajr(c — o) r(a)r(o) 



Invoking (31 ), 6 > 0, and a + 6 — c + 1 = r^br, we deduce that 
P*(r)/x(r) ~ 2 r(c)r i° + ^" C+1) (l-r 2 )^ + --^ ~ 1 as < 1 - r « 1, 

r(o)r(6) 



and thus (22) is violated. Consequently, f\ is not in V(7ii), an d thus not an eigen- 
function of Tie. 

As the above linear transformation has a pole when a + b — c £ *Z, we have to use 
substitute formulas to investigate the limiting behavior at z — 1. Such transforma- 
tion formulas can be found in [U 15.3.10, 15.3.12]. Instead of displaying the explicit 
expressions in the sequel, we will just quote the limiting behavior of F(a, b; c; z) as 
z converges to 1. We note that c < a + b+1 since a + b — c+1 = —^—r > 0. We have 

o m— 1 

F(a, b; a + b; z) « --fr^TTT ln(l - 2) as < 1 - z < 1, 

r(o)r(6) 



and for fc G N, 



r^rco + ft-*),, ,_ fc 



F(a, b; a + b - k; z) « ^-7 — — -r -(1-^) as < 1 - z < 1. 

r(o)r(6) 

For every fcsNU {0}, the case a + b — c = k is equivalent to -— r — k + 1, and thus, 
a short computation using (31) yields 

P*(r)f[(r) ~ (1 - r 2 )^r- (fe+1) = 1 as < 1 - r < 1. 



This violates (22) and thus f\ is not an eigenfunction of Hi. 

Step 3. For —a G N U {0} ; it holds f\ G T>{^H.(), and thus f\ is an eigenfunction of 
Tii. The corresponding eigenvalue is X = £ + Ik + Ik (k + £ + y — l) (m — 1). 
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We set a = —k for some fceNU {0} and observe that the hypergeometric series in 



(24) terminates, namely 



It follows that /i is a polynomial of degree 2k + £ and thus a continuous function 



on [0, 1]. In particular, the asymptotic boundary condition (22) is trivially satisfied. 
Regarding the integrability of f\, we only have study the singularity if £ > 1. 
In this case, we have that {f\(r)) 2 r N ~ 3 « r 2^+JV-3^ -^j^ j s integrable, and thus 
/i G i?£. Moreover, since Hefi — A/i we also know that 7^/i G i7|, and thus 
/i G V(Hi)- We conclude that /i is an eigenfunction of %£ if a is a nonpositive 
integer. 

It remains to express the eigenvalue A in terms the fixed variables. A straight forward 



computation based on the expression (27) yields that 



X = £ + 2k + 2k(k + £ + 1 J (m-1). 

This concludes the proof of Proposition [3] □ 

Appendix: Elliptic theory in Hi 

This appendix provides some results related to elliptic problems in Hi that apply 



in the discussion of the displacement Hessian 7i in Section 2J2_. The statements in 
the Lemmas [2]-[5] below are the analogs of very classical facts in standard Sobolev 
theory: density of smooth functions, a Hardy-Poincare inequality, a Rellich lemma, 
the Poisson problem. Our Hilbert space i?* differs from the well-studied Sobolev 
space H x {Bi) only by the weight p*, which is, by the regularity of the Barenblatt 
profile in the ball B\, a very mild variation: p* is finite and degenerates only on the 
boundary dB\. In this regard, it is not surprising that the results in H^ can be 
proved analogously to or based on the known theory. 

We first give an overview of the main facts. We start with a classical density state- 
ment. 

Lemma 2 (Density of smooth functions). For every ip e H 1 , there exists a sequence 
{"0v}i, e N ^ n C°°(Bi) such that 



lim / p*|VO-Vv)|Gte = 0. 

!/f00 J 

The next result is an embedding theorem between weighted Lebesgue and Sobolev 
spaces, cf. [HI Corollary 13]. 
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Lemma 3 (Hardy-Poincare inequality). Let p > 1 such that p + m > 3. Then 

in^j ff- m ^-c) 2 dx < j p*\ViP\ 2 dx (34) 

for any ip G i?* . 

Notice that the infimum in the statement of the Hardy-Poincare inequality is at- 
tained by c = J pP~ m ip dxj J p v ~ m dx. The inequality is sharp if p + m — 3. In most 



cases, we will apply (34) with p = 2. The general result is used in the proof of the 
following compactness result. 

Lemma 4 (Rellich lemma). The embedding of H l in L 2 {Bi,p 2 ~ m dx) is compact. 

We finally study the Poisson problem in H* . 

Lemma 5 (Poisson problem). For every u G L 2 (Bi, p™~ 2 dx) with J udx = there 
exists a unique (up to an additive constant) ip G H 2 oc (Bi) fl H x such that 

— V • {p*Vip) — u in Bi 
p*Vip ■ v — on dB\. 

As before, we understand the boundary conditions in the above lemma in the sense 



of (17). 



We first address the proof of Lemma |2j which uses classical approximation tech- 
niques, see e.g., [HI Sec. 5.3]. 

Proof of Lemma^ We fix ip G H 1 and notice that that ip G H X {K) for every 
compact set K C B\ since p* is positive and finite away from dB\. In a first step, 
we show that ip can be approximated by C°°(Bi) functions. For this purpose, we 
consider An = B-i and annuli A v = B»+2 \ B_^_ for every v G N, where B r denotes 

3 i/+3 "+1 

the ball of radius r around the origin, to the effect of Bi = [J^ =0 A v . Let {^}„ e]No 
be a partition of unity subordinate the covering {A u } ue¥So and let {C} mG n be a 
sequence of standard mollifiers. We fix e > arbitrarily. Since r] v ip G H 1 (A U ), for 
any v G N there exists a /i = p(e, z/) G N such that 

\ 1/2 

P*\V{r] u ip - (r) v ip) *([+) \ 2 dx) < 



2 u+l 

and 

supp ((^V) * O c ^^ \ -Bh=i- 



We define ip = Yl^Lo^^ where ip v = {r] v ip) * (^. We obviously have ip G C°°(B 1 ] 
Moreover, since ip = Yl^LoVf^^ we have for every compact set K C B\\ 



1/2 oo , r \ 1/2 



f P *\V(iP-iP)\ 2 dx) < ^(jp,\V{r] v iP-iP v )\ 2 dx\ 



v=0 

oo 
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Taking the supremum over all K, this shows that C°°(Bi) is dense in H*. 

In order to prove density of functions that are smooth up to the boundary, we may, 
due to the above argumentation, suppose that ip G C°°(B\). For x G B\ and i/eN, 
we define ip v {x) = ip((l — \)x). It follows that ip v G C 00 (B 1 ) and 

p*\V{Tp-i> v )\ 2 dx U -H 

by the dominated convergence theorem. □ 

We now turn to the proof of Lemma [3j which relies on the following one-dimensional 
inequality: 

Lemma 6 (Double- weight Hardy inequality). Let p > 1 such thatp + m > 3. Then 

inf / {\-r)^{Tl)-c) 2 r N - x dr < f (1 - r)^i(^) 2 ^" 1 dr (35) 

c6R Jo Jo 

for any ip G C^O, 1]. 

Notice that the above inequality is optimal for p + m = 3 in the case m/2. For 
p — 1 and m — 2, however, this estimate fails (logarithmically) as the underlying 
Hardy inequality is critical. 

Proof of Lemma^ Since C°°{B{) is dense in B 1 by Lemma ([l|, it is enough to 
consider the case where ip is smooth up to the boundary. We prove a slightly 
stronger statement by choosing c = f pP~ 1 \x\~ 1 ipdx/ f pP^ 1 ^^ 1 dx, or equivalently 
by assuming that 

- x \x\~ l i)dx = 0. (36) 



Notice that the integral is well-defined for every ip G C°°(Bi) by the assumption on p 
and since N > 2. We rewrite the statement in spherical coordinates: Let x = ru and 
ip(r,u) = ip(x) where r = \x\ and uj = x/\x\. Then |\7-0| 2 = {d r ip) 2 + -^\V$n-iiIj\ 2 , 



where V$jv-i denotes the tangential gradient and thus, by the definition of /?*, (34) 
follows from 

I \ P~ m ^O AT 1 

/ (l—r) m - 1 i/jr ~ dudr 
il Js N ~ 1 

< if (1 - r)^i ( (d r ij) 2 + \\ VsN-iijA r N ^ dudr. (37) 



r*J 



o J$ 



N-i \ r 



For any fixed u>, we consider ip^^r) = ip(r,u). Since ip u G C°°[0, 1], we can apply 
Lemma M componentwise for any to G § N ~ 1 . Observe that the infimum in (35) is 



attained by a constant c ~ J (1 — r) '"- 1 if) u r N x dr, so that (35 ) and integration over 
S N ~ l yield 
i 

^ ! ./n 



ls N - 1 \Jo J 7$jv-i Jo 



(l-r)^=iipr"- L dr ) du + / / (1 - r)^(d r V>) r drdw. 
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To derive (37), we need to control the first term on the right. Integrating by parts 



and applying the Cauchy-Schwarz inequality, we compute 



(1 — r)™- lr ipr dr 

"'- 1 (N-1) f\\-r)^ijr N - 2 dr+ 1 ^^ f\l-r)^d r i>r N - x dr 



< 



P~ 1 JO 

i 



p — 1 



2-1 



(1 — r)™- 1 ^ dr 



ii 



rl 2 3 \ 1/2 / ,i 

( / (l-r^r^dr 



(l-r)^^) r dr 



1/2 



Now that the above estimate becomes 
-l 

2„JV-1 



s^- 1 Jo 



< 



(l-rj^^V-Mrrfw 



(1 — r) m -i , 0r 2 dr ) du; + 



s^- 1 Jo 



(l-rj^^r" 1 *^ 



and the statement in (37) follows from the estimate 



$N-1 \J0 



(l-r)^ipr N - 2 dr) du < 



s^- 1 Jo 



1 - r)^\V S N-iilj\ 2 r N - s drdu. 

(38) 

To prove (38), we start by recalling the Poincare inequality on a sphere (cf. [14, 
Theorem 2.10]) 



edu < 



£duj) + / \V$N-i£\ 2 du, 



P-i 



which we apply to C, = L (1 — r) m -! ^r^ 2 dr to the effect of 



i p _ i \ 2 

(1 — r) m ~! %jjr N ~ 2 dr J dw 
s- v ! V./u / 

1 \ 2 

(l_ r )^^i^ r ~ drdu \ 

■■: ■ / I / (1 -r)^|V S iv-i^|r JV " 2 dr ) dw. 

s*- 1 Vio / 



< 



The first term on the right vanishes since 



s^- 1 Jo 



(l-r)^itjjr N - 2 drduj ~ p^lx^dx ^ 



(>!(,) 



2N 



The second term can be estimated using the Cauchy-Schwarz inequality in the inner 
integral: 



p-i 



1 -r)^|V S iv-n/V ^ r ) duj 

1 2 _ 3 /■ /-I 

< / (1 — r) : ^ I r N ~ 1 dr 



o 



s^- 1 Jo 



1 - r )^ r r|V$iv-i^|V JV " d rfro?Cc;. 



Since the prefactor is bounded by the assumption on p, this proves (38). 



□ 



Proof of Lemma^ The statement (35) follows as a combination of the two Hardy 
inequalities 

•1/2 r l/2 

(V>-V(V 2 ))V v ~ 1 dr < / (V^V^dr 



and 



Indeed, 



i /-l 



(l-r)'^(V'- -0(1/2)) or < / (l-r)^T(V>) dr. 

1/2 A/2 



(39) 



(40) 



(1 _ r )^=i{^ - ijj{l/2)Yr N -'dr 



< 



1/2 pi 

(if; - if;(l/2)) 2 r N ~ l dr + / (1 - r)^(i/j - 0(l/2)) 2 dr 



o 



1/2 



p9l&p0l /-1/2 



< 



li A/2 

1 

i\2„N-l 



(l-r)^(VO r *■, 



o 



and the term on the left can be bounded from below by taking the infimum over all 
c = -0(1/2) eR. 



For proving (39) and (40), it is enough to consider the case where "0(1/2) = 0. The 



proof of the Hardy inequality (39) is standard: 

1/2 1 ._ r=1 / 2 2 f 1/2 



if?r"- 1 dr 



N 



tfj 2 r N 



r =o *Jo 



ipip r dr 



< If [ 1/2 ^r"-idr 



N 



1/2 







1/2 



1/2 



(V/) V v+i rfr 



We apply the Young inequality 2a6 < a + b to infer ( 39 ) . 
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Likewise, for (40) we have 

• i 



p — m 



(1 — r)™- 1 ?/) dr 

1/2 



m — 1 P-1 r 

(1 — r) m - 1 z/>'' 

p- 1 V ; v 



r=\ 1 »l 

m — 1 / . 



+ 2 / (1 — r) rn - 1 ipip' dr 

r=l/2 P " X A/2 



< 2— I {l-r)^ip 2 dr\ ( (1 - r)^ (if)') 2 dr 



■>.,•■ ■■', 



In order to deduce (40), it remains to check that (1 — r)™- 1 < (1 — r)™- 1 holds if 



p > 3 — m, and to apply the Young inequality. □ 

We now address Lemma HI 

Proof of Lemma [7} We derive the statement from the well-known analogous state- 
ment for regular Sobolev spaces, cf. [UJ p. 272]. Let {^„} ne N denote a bounded 
sequence in H*. We may without lost of generality assume that J p^~ m ip n dx = 0, 
where p < 2 denotes a constant that satisfies the hypothesis of Lemma [3] (this is pos- 
sible since m > 1), and thus the sequence is uniformly bounded in L 2 {B\^p v ~ m dx) 



by (34). As H l pt is a separable Hilbert space, we may extract a subsequence that 
converges weakly to a limit if) both in H x and in L 2 {B\^p v ~ m dx). Moreover, as p* 
is bounded below by a positive constant in every compact subset of B\, we have 
that the sequence {V'nlneN is bounded in the (unweighted) Sobolev space H 1 (B 1 _i) 

k 

for every k e N. By the Rellich compactness lemma, we may successively extract 
further subsequences {VVj,(n)}neN that converge to if) strongly in L 2 (B l _i). By the 
boundedness of p*, we may choose a diagonal sequence r(k) = Tk(N(k)) such that 

pt- m W-A(k)) 2 dx < I, (41) 

B, i K 

1 fc 

for all k e N. For e > arbitrary but fixed, we write 

P 2 - m ^-^ T{k) fdx = ! P i-™^-^ T{k) y d x+ ! ^ P 2 - m {^-^ T{k) ) 2 dx. 



The first integral converges to zero if k goes to infinity thanks to (41). For the 
second one, we have that 



P l~ m ^-^ { k)) 2 dx < £^T / ffty-^fdx, 



2-P 

V(fc)J ax - fc™- 1 / p; ^-^T(fc) 

Bi\Bi_ e ^Bi 

and integrals on the right are bounded uniformly in k. Hence, letting k converge to 
infinity, we have that 



lim [ pl- m ^-^( k )) 2 dx < Ce 



m—1 

) 



for some C > 0. Since e was arbitrary and p < 2, this proves the statement of 
Lemma HI □ 
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The proof of Lemma [5] is very classical and we just sketch it. 

Sketch of the proof of LemmaUn Observe that the elliptic problem is the Euler- 
Lagrange equation of the strictly convex functional 

•^""(VO = o / P*|V?/>| 2 <ix — / uip dx 

for ip G if* . Therefore, it suffices to prove existence of a minimizer of J 7 in if * . In- 
deed, by the strict convexity of J 7 , minimizers if) are unique up to additive constants 
and satisfy the weak Euler-Lagrange equation 



/ p*V^ • V^dx = I u£ dx for all f G H l pf . 



(42) 



In particular, choosing £ G C%°(Bi), we see that ip is a distributional solution of 
— V-(p*V^>) = u, which also holds in the strong sense by interior regularity estimates. 



Moreover, with this information, (42) becomes 



J p.V^p-Vidx = - ! V • (p* WO f rfx for all f G #^ , 

i.e., p^Vip ■ v = on <9f?i. 

To prove existence of minimizers of T in H \ t , we just mention the basic steps, 
following the direct method. We consider a minimizing sequence for J 7 in H^. 
With the help of Lemma pi it is easy to check that this sequence is bounded in H^ t . 
This fact is based on the estimate 



h dx = /'""-^ 



/ r \ 1/2 / r \ 1/2 

< [J p?- 2 u 2 dx\ U fir m ty-c) 2 dx 

< (Jp r r 2 u 2 dx^ (| p ,iv^i 2 ^ , 

where c denotes the optimal constant of Lemma [3j that could be introduced in the 
first identity thanks to the fact that u has average zero. As a consequence, we can 
find a weakly converging subsequence. By lower semicontinuity of iF with respect 
to weak convergence, we immediately deduce that the minimum of iF in if* is 
attained. □ 
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